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THE  REBATRON  AS  A  HIGH  ENERGY  ACCELERATOR 


I.  Introduction 

The  electric  field,  that  is  responsible  for  the  acceleration 
of  charged  particles  can  be  written  as  the  sum  of  an 
electrostatic  and  an  inductive  component.  The  electrostatic 
component  is  generated  by  electric  charges  and  its  divergence  is 
proportional  to  the  charge  density  while  its  curl  is  zero.  The 
Inductive  electric  field,  on  the  other  hand,  is  generated  by  a 
time  varying  magnetic  field  and  its  curl  is  proportional  to  the 
rate  of  change  of  magnetic  field,  while  its  divergence  is  zero. 

It  is  the  inductive  field  that  is  responsible  for  the 
acceleration  of  particle  In  Induction  accelerators^”^. 

Only  modest  acceleration  can  be  achieved  with  electrostatic 
fields  because  the  maximum  potential  is  limited  by  insulator 
flash-over.  However,  this  is  not  the  case  with  inductive  fields, 
since  the  voltage  is  induced  only  in  the  circuit  threading  the 
flux  and  the  voltages  of  a  sequence  of  circuits  successively  add 
to  the  energy  W  of  the  accelerated  particle.  The  striking 
difference  between  the  two  fields  becomes  apparent  when  the  orbit 
of  the  particle  is  circular.  After  a  complete  revolution  W2  -Wj 
a  0  for  an  electrostatic  field,  but  Wg  -  *  0  for  an  inductive 

field,  if  the  magnetic  flux  through  the  particle  orbit  changes 
with  time.  An  additional  advantage  of  induction  accelerators  is 
their  inherent  low  impedance.  As  a  result,  these  devices  are 
ideally  suited  to  drive  high  current  beams. 

Quite  naturally.  Induction  accelerators  are  divided  into 
linear  and  cyclic.  The  linear  devices  are  in  turn  divided  in 
Astron-type®”^,  Radi ac-type^^”^^ ,  and  auto-accelerators 

Manuacript  approved  November  27, 1986. 


1 


In  the  first  type,  ferromagneti c  induction  cores  are  used  to 
generate  the  accelerating  field,  while  "air  core"  cavities  are 
used  in  the  second.  In  the  auto-accelerator  the  air  core 
cavities  are  excited  by  the  beam's  self  fields  rather  than 
external  fields.  Similarly,  cyclic  accelerators  can  be  divided 
into  two  categories.  In  the  first  category  belong  those  cyclic 
accelerators  that  use  weak  focusing  to  confine  the  electron  beam 
as  in  the  con ventional and  the  modified  betatron^^"^®.  In 
the  second  category  beam  confinement  is  achieved  with  strong 
focusing  fields.  This  latter  category  includes  the  Stellatron*  , 
Racetrack^^,  and  Rebatron  accelerators^^.  All  three  devices  are 
currently  in  a  rather  preliminary  state  of  development. 

In  the  last  two  devices^^*^^  the  particle  acceleration  occurs 
rapidly,  i.e.,  in  a  few  usec  and  thus  limitations  imposed  by 
instabilities,  field  errors  and  synchotron  radiation  losses  are 
substantially  relaxed. 

2  1 

A  rebatron  is  shown  schematically  in  Fig.  1.  The  high 
gradient  localized  field  that  is  responsible  for  the  rapid 
acceleration  is  produced  by  a  convoluted  parallel  transmission 
line,  although  other  transmission  lines  may  be  more  appropriate 
in  an  actual  system.  Since  the  acceleration  occurs  over  a 
few  usecs,  the  contraints  imposed  on  the  vertical  field  are  very 
stringent.  To  reduce  the  inductance  of  the  system,  the  vertical 
field  is  generated  by  two  coaxial,  cylindrical  plates  as  shown  in 
the  lower  right  corner  of  Fig.  1.  The  axes  of  these  plates 
coincide  with  the  major  axis  of  the  toroidal  vessel  and  are 
located  symmetrically  around  the  minor  axis  of  the  torus.  The 
purpose  of  the  gap  in  the  outer  plate  is  to  provide  a  field  with 
the  desirable  external  field  index.  The  transmission  lines 
change  mainly  the  local,  vertical  magnetic  field,  while  the 
magnetic  flux  through  the  beam  orbit  remains  approximately 
constant . 


The  mismatch  between  the  beam  energy  and  the  vertical  field 
is  alleviated  by  a  strong  focusing  field.  This  field  is 


generated  by  a  set  of  t  *  2  torsatron  windings,  i.e.,  two  twisted 

wires  that  carry  current  in  the  same  direction.  In  addition  to 

the  transverse  components  of  the  field,  the  torsatron  windings 

provide  a  zero  order  toroidal  magnetic  field.  The  beneficial 

effect  of  twisted  quadrupole  fields  on  the  beam  orbits  has  been 

2  2 

recognized  for  several  years  .  In  recent  years,  the  t  *  2 

stellarator  field  has  been  used  by  Roberson  et  al.  to  Improve  the 

bandwidth  of  the  modified  betatron  in  a  device  named  Stellatron^^ 

20 

and  also  in  the  racetrack  accelerator  ,  a  device  similar  to 
rebat  ron . 

2 1 

In  recent  studies  ,  we  have  demonstrated  that  the  energy 
bandwidth  of  a  rebatron  accelerator  can  be  very  wide.  In 
addition,  we  have  shown  that  the  torsatron  windings  also 
substantially  improve  the  current  carrying  capabilities  of  the 
device  and  could  alleviate  the  beam  displacement  problem 
associated  with  the  diffusion  of  the  self  magnetic  field. 
Therefore,  it  is  expected  that  the  limiting  current  in  a  rebatron 
would  be  determined  from  collective  instabilities  and  not  from 
the  macroscopic  stability  on  the  beam  orbits. 

Beam  capture  in  the  rebatron,  as  in  other  devices  that  use 

strong  focusing  fields,  is  very  difficult.  The  reason  is  that 

the  strong  focusing  fields  make  the  particle  orbits  insensitive 

to  the  energy  mismatch  and  thus,  small  changes  in  the  vertical 

magnetic  field  are  not  sufficient  to  move  the  beam  from  the 

injection  position  near  the  wall  to  the  minor  axis  of  the 

torus.  Recently,  we  have  developed  two  injection  schemes  that 

appear  very  promising.  One  is  based  on  the  drag  force  of  a 

2  3  2  4 

resistive  chamber  wall  and  the  other  on  the  modification  of 
the  beam  orbit  by  a  time  varying  field. 

During  acceleration,  when  the  vertical  magnetic  field 
exceeds  by  far  the  torsatron  field  B^,  the  beam  dynamics  is 
solely  determined  by  the  vertical  and  toroidal  fields.  For  most 
applirations  of  interest  B^  >>  B^  at  the  peak  of  acceleration  and 


thus,  beam  extraction  from  the  rebatron  is  similar  to  that  from  a 

9  1 

modified  betatron^^. 

In  this  report,  we  present  results  from  our  studies  of  the 
beam  acceleration  in  a  rebatron  accelerator  when  the  magnetic 
fields  varies  with  time.  In  section  II  the  applied  electric  and 
magnetic  fields  are  described  and  in  section  III  a  simple 
theoretical  model  is  examined  and  some  important  conclusions  are 
derived  from  this  model.  The  numerical  results  for  an 
accelerated  beam  and  their  interpretation  are  given  in  section 
IV.  Finally,  the  conclusions  of  this  investigation  are  presented 
in  section  V. 


1 1 .  Appl led  Fields 

a.  Torsatron  and  External 
Toroidal  Magnetic  Fields 

In  the  local  cylindrical  coordinate  system  e  ,  e.,  e 

P  <9 

shown  in  Fig.  2,  the  magnetic  field  components  of  the 
torsatron  and  external  toroidal  coils  are  given  by  the 
expressi ons  ^ 
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The  coefficients  and  A^^^  are  equal  to 
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The  remaining  parameters  are  defined  as  follows: 
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Here,  is  the  current  flowing  in  the  windings,  is  the 

winding  radius  and  L  is  the  period  of  the  windings,  2  ^  is  the 

width  of  the  current  carrying  conductor,  r  is  the  major  radius 

ext  ^ 

of  the  torus,  is  the  external  toroidal  magnetic  field  at  the 

minor  axis  of  the  torus  and  ^^(x), 

modified  Bessel  functions  of  the  first  and  second  kind  and  their 
derivatives.  In  a  toroidal  device,  the  period  satisfies  the 
rel ati on 


where  N  is  the  number  of  turns  around  the  torus  for  each  winding. 

The  zero  order  fields  B^®^,  b[®^,  are  the  field  components 

p  <0  s 

produced  by  the  helical  windings  in  a  straight  (cylindrical) 
configuration  while  the  terms  proportional  to  o r ^  are  the  first 
order  toroidal  corrections.  These  fields,  as  given  by  Eqs.  (2.2) 
and  (2.3),  have  been  obtained  for  the  surface  current  density 


J  *  0 

13 


(2.10a  ) 


I  aOg 

2 ~I  +  (o/r^)  cos 4  “S), 


(2.10b) 


25o„ 
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where  f(4)  is  a  periodic  function  of  6  with  period  t,  and 


lif-5<4<5 
^^'*'^*0if  6<4<it-5 


(2.11) 


From  Eqs.  (2.3d)  and  (2.3e)  it  follows  that  there  is  a 
constant  component  of  the  torsatron  magnetic  field.  This 

component  is  due  to  the  fact  that  the  currents  in  the  two 
torsatron  windings  are  in  the  same  direction.  (In  a  stelarator 
this  component  is  zero).  It  is  assumed,  in  the  rest  of  this 
paper,  that  this  undesirable  component  is  cancelled  by  adding  an 
external  magnetic  field. 

The  three  magnetic  components  at  s=0  are  plotted  in  Figs.  3a 
and  3b,  for  4  =  0  and  4  *  ir/2,  respectively.  The  various 
parameters  for  these  plots  are  listed  in  Table  I  and  are  those 
used  for  the  torsatron  and  external  toroidal  magnetic  field  in 
the  following  sections.  At  4=0,  the  8^  component  is  equal  to 
zero.  The  3^  component  increases  linearly  with  p  near  the  minor 

<P 

axis  and  much  faster  near  the  windings.  The  toroidal  correction 
at  p  =  0  is  approximately  -100  Gauss.  In  Fig.  3,  in  addition  to 
the  toroidal  field  produced  by  the  torsatron  windings,  there  is 

P  y  ^ 

an  external  toroidal  field  equal  to  10  kG  at  p  =  0  ,  which 

is  produced  by  a  set  of  toroidal  coils.  As  seen  from  Eqs. 


I 


(2.1c)  and  (2.2c)  this  field  varies  as  1/r  and  Is  Independent  of 
the  toroidal  angle  e  *  -  s/r  .  At  p  =  0  ,  the  total  toroidal 

a  y  ^  ^ 

field  BgQ  +  B^*  15.6  kG.  In  general,  B^  does  not  vary  as  1/r, 
since  it  depends  on  the  toroidal  angle. 

Figure  3c  gives  the  magnetic  field  lines  in  the  (r,z)  and 
(p,s)  planes.  Because  of  the  toroidal  corrections,  the  magnetic 
axis  does  not  coincide  with  the  minor  axis  of  the  torus,  located 
at  r=100  cm.  The  magnetic  axis  is  always  shifted  toward  the 
major  axis  of  the  torus.  This  shift  is  due  to  the  contribution 
in  the  magnetic  field,  at  some  point  inside  the  torus,  mainly 
from  that  section  of  the  windings  which  is  located  at  a  toroidal 
angle  180“  from  this  point. 

In  the  theoretical  model  that  is  formulated  in  Section  III, 
the  toroidal  corrections  are  neglected  as  well  as  all  the  terms 
with  m  >  2  in  Eq.  (2.1)  through  (2.4).  Also,  it  is  assumed  that 
-5  <<  1,  so  that  C,  =«  2.  Finally,  use  is  made  of  the  cylindrical 

•  AAA 

components  B^,  B^,  B,  in  the  global  coordinate  system  e^,  e.,  e_ 
shown  in  Fig.  2.  Their  relation  to  the  magnetic  components  in 

AAA 

the  local  coordinate  system  e^,  e^,  e^  is  as  follows 


B  cosi)  -  B  .  sin*, 
0  ^ 


8  sin*  +B^  cos*, 
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Be  =  -  8s 


(2.12a) 
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Under  these  conditions,  the  magnetic  field  from  the  torsatron  and 
the  toroidal  coils  is  approximately  equal  to 
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where 
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^1  “  Oq  ><2  (og),  (2.14a) 

Eqs.  (2.1)  through  (2.8)  were  used  to  compute  the  magnetic  field 
in  the  numerical  integration  of  the  equations  of  motion  for  the 
beam  centroid.  In  this  case  too,  it  was  assumed  that  ^  <<  1,  so 
that  -  2  and,  in  order  to  save  computer  time,  only  one  term, 
namely  m  *  1,  was  retained  in  the  series. 

b.  Vector  Potential  and  Magnetic  Field  of  the  Plates. 

Let  a  single  cylindrical  plate  of  width  b  and  height  h, 
where  b  <<  h,  be  located  so  that  the  axis  of  the  cylinder 
coincides  with  the  major  axis  of  the  torus  and  the  distance  of 
the  plate  from  the  minor  axis  of  the  torus  be  equal  to  a.  For  a 
homogeneous  current  density  with  a  single  component 

=»  Ip/(bh),  where  Ip  is  the  current  flowing  in  the  plate,  the 
vector  potential  is 


where 
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The  point  (x,z),  where  A  is  computed,  is  In  the  local  coordinate 
system  e^,  e^,  e^  of  Fig.  2,  and  K'(k'),  E'(k")  are  the 
complete  elliptic  Integrals.  When  the  plate  Is  located  close  to 
the  torus,  so  that  a/r^  <<  1,  and  one  Is  Interested  In  the  vector 


potential  only  inside  the  torus,  then  ,,^2..  ,  and  the  elliptic 

(26)  2 

integrals  can  be  expanded'  '  with  respect  to  k'  .  To  lowest 

A 

order  A^  is  given  by 
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which  contains  the  toroidal  corrections  to  the  lowest  order  in 

a/r^.  A  further  si mpl i cati on  occurs  from  the  fact 

that  b  <<  h,  so  that  A.  is  computed  in  the  limit  b  +  0.  Then  the 

integrals  in  Eq.  (2.17)  can  be  computed  exactly  and  A^  is  equal 

to 
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The  Bp,  Bj,  components  of  the  magnetic  field  near  the  plate  (and 
therefore  inside  the  torus),  and  to  the  lowest  order  in  a/rg,  are 
equal  to 
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where 
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Substitution  of  Eq.  (2.18)  into  Eqs.  (2.20a),  (2.20b),  leads  to 
the  expressions 


h  2  2 

1  la-y  (j+Z)  +  (a-x) 

F  - 5 - r 


2.21a) 


(f  -  z)  +  (a  -  x) 


and 
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In  the  limit  r^  +  •,  Eqs.  (2.20),  (2.21)  reduce  to  the  magnetic 
field  of  an  infinite  length  slab,  as  they  should. 

The  particle  orbits  are  not  stable  unless  the  external  field 
index  n  of  the  vertical  magnetic  field  is  in  the  region 
0  <  n  <  1.  Such  a  field  can  be  generated  by  a  set  of  two  plates 
which  surround  the  torus  (cf.  Fig.  1),  with  the  outer  plate 
having  a  gap  symmetrically  located  in  the  middle  of  the  plate. 
From  the  vector  potential  for  a  single  plate,  it  is  easy  to 
compute  the  vector  potential  of  the  outer  plate  with  a  gap  of 
height  H^.  This  vector  potential  is  given  by 


h  -  H  a, h,0)  _  l_l  A^(x,z,r^,a,H^,0).  (2.22a) 


0’  ’  0 


The  current  density  in  the  presence  of  a  gap  is  proportional 
to  In/(h  -  H  ).  The  first  term  in  Eq.  (2.22a)  is  the 
contribution  to  from  the  whole  plate  (of  height  h)  and  the 
second  term  is  the  contribution  of  the  gap  (of  height  H^)  which 
should  be  subtracted.  Since  the  inner  plate  is  located  at  a 
distance  -a  from  the  minor  axis  and  carries  a  current  -  Ip,  its 
vector  potential  is 


Ag(x,z,rg,-a,h.O)  , 


and  the  combined  potential  from  both  plates  is 


(2.22b) 


(2.23) 


The  magnetic  field  is  computed  with  the  help  of  Eqs.  (2.13) 
through  (2.23).  Near  the  minor  axis,  the  magnetic  field 
simplifies  to  the  following  linearized  expressions: 
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Equations  (2.24a),  (2.24b)  are  Identical  in  form  to  the 
linearized  betatron  field  in  the  modified  betatron.  According  to 
Eq.  (2.25a)  the  external  field  index  at  x®z=0  is  sensitive  to 
the  height  of  the  gap.  In  particular,  if  there  is  no  gap  on 
the  outer  plate  (i.e.,  H  »  0),  then  n  becomes  negative. 


The  magnetic  field  components  of  the  plates  are  plotted 
without  a  gap  in  Fig.  4  and  with  a  gap  on  the  outer  plate  in  Fig. 
5.  The  various  parameters  for  these  plots  are  listed  in  Table 
II.  The  effect  of  the  gap  on  the  magnetic  field  is  quite 
prominent.  In  the  absence  of  a  gap,  n^  changes  from  positive  to 
negative  as  r  increases  from  90  to  110  cm  for  z  less  than  10 
cm.  On  the  other  hand,  when  the  gap  is  present,  n^  has  the 
desired  value  of  0.5  near  the  minor  axis  of  the  torus  (i.e.,  in 
the  linear  regime),  but  as  one  approaches  the  gap,  n^  varies  from 
positive  to  negative  as  z  increases  from  zero  to  10  cm.  Figures 
4b,  5b  and  4d,  5d  show  the  variation  of  B^  versus  r  and  B^  versus 
z,  respectively,  for  various  values  of  z  or  r.  Ideally,  they 
should  be  straight  lines,  parallel  to  the  horizontal  axis  in  the 
region  of  interest.  As  may  be  seen  from  Fig.  5  this  is  the  case 
when  there  is  a  gap.  Although  only  the  first  order  toroidal 


corrections  were  retained,  they  are  quite  noticeable.  For 
example,  in  Fig.  4a  at  z^O,  there  is  a  difference  of  0.6  Gauss/kA 


in  the  values  of  at  r=90  cm  and  r*110  cm,  due  to  the  toroidal 


correction.  Finally,  it  should  be  pointed  out  that,  because  of 
the  gap,  the  value  of  the  vertical  magnetic  field  at  the  minor 
axis  drops  from  18.45  Gauss/kA  to  17.58  Gauss/kA. 


As  already  mentioned,  in  order  to  keep  the  beam  matched,  the 
current  in  the  plates  should  Increase  in  synchronism  with  the 
energy  of  the  beam.  Therefore,  there  is  an  electric  field 


Ep  associated  with  the  time  dependent  vector  potential  of  the 
plates,  namely 


1_ 

c 


3  A 


3  t 


(2.26) 


As  shown  in  Section  III,  this  electric  field  will  have  an  effect 
on  the  rate  of  change  of  the  current  in  the  plates. 


C.  The  Accelerating  Electric  Field. 

Consider  a  periodic  sequence  of  cylindrical  perfectly 
conducting  sections  of  radius  with  their  axes  lying  along  the 
same  line  and  separated  by  a  distance  d,  as  shown  in  Fig.  6a. 

Each  alternate  section  is  charged  to  a  potential  +  V^/2  and  the 
remaining  sections  are  charged  to  a  potential  - 
period  of  the  structure  is  equal  to  0,  where  D  >>  d.  Under  this 
assumption,  the  local  electric  field  in  the  region 
-  D/4  <  s  <  0/4  is,  for  all  practical  purposes,  the  same  as  that 
of  two  cylinders  of  infinite  extent.  In  either  case,  the  local 


electric  field  is  given  by  the  solution  of  Laplace's  equation, 
i.e.,  A  V  =  0.  The  advantage  of  solving  this  equation  for  the 
periodic  structure  rather  than  the  two  cylinders  of  infinite 
extent  is  that  it  leads  to  a  solution  expressed  in  terms  of  a 
much  faster  converging  series.  The  lowest  order  toroidal 
corrections  have  been  included  in  the  solution.  The  method  for 
incorporating  these  corrections  is  the  same  as  that  used  for  the 
torsatron  magnetic  field.  The  operator  a  in  the  local  coordinate 
system  e  ,  e  ,  e.  can  be  written  as  follows: 


A  = 


+  5A 


(2.27) 
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and  s  *  -  r^e.  The  potential  V  is  also  written  as  a  sum  of  the 
zero  order  and  first  order  contribution  in  the  toroidal 
corrections,  namely,  V  *  where  satisfies  the 

zero  order  Laplace  equation 


a(»>  v*"*  .  0 


(2.29) 


Since  the  orginal  Laplace  equation  becomes  equal  to 


(A'""'  +  5A)  =  0. 


(2.30a  ) 


it  follows  from  Eq.  (2.29)  that  the  first  order  toroidal 
correct  ion  in  the  potential  satisfies  the  equation 


i(»)  yd)  =  -  5.  y(»). 


(2.30b) 


The  method  of  computing  the  toroidal  corrections  consists  in 
solving  first  Eq.  (2.29)  with  the  boundary  condition  at  p  = 


0  <  s  <  j 
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(y-  -  s )  if 
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and  V^^^(-s)  »  -  (s).  Then  is  computed  by  solving  Eq. 

(2.30)  with  the  boundary  condition  that  at 
0  »  p^ ,  0,  if  -  D/2  <  s  <  D/2.  The  accelerating 

electric  field,  being  the  gradient  of  -V,  is  computed  from  the 
rel ati ons 
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The  components  of 

the  accelerating  electric  field. 

including 

the  lowest  order  toroidal  corrections,  are 

0  0  0 

(2.33a) 
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The  rest  of  the  parameters  in  these  relations  are  defined  as 
foil ows : 
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The  electric  field  is  computed  in  the  interval  -D/4  <  s  <  D/4. 

A  choice  of  D/4  •  3p^  leads  to  a  fairly  accurate  representation 
of  the  accelerating  field.  As  a  result  of  the  modified  Bessel 
functions,  the  series  in  Eqs.  (2.34)  and  (2.35)  converge 
considerably  faster  than  those  used  previ ously ^ .  It  may  also 
be  seen  from  these  equations  that  as  o  increases  from  zero 
to  p^  more  terms  must  be  retained  in  the  series  for 
convergence.  Finally,  the  presence  of  the  toroidal  corrections 


and,  in  particular,  the  1/r  dependence  of  the  s-component  of  the 
electric  field  (cf.  Eq.  (2.33c)),  is  essential  in  the  orbit 
equations.  It  will  be  shown  in  Section  III  that  the  toroidal 
corrections  of  the  accelerating  electric  field  cancel  out  the 
first  order  contribution  from  the  electric  field  of  the  plates 
(Eq.  (2.26))  in  the  expression  for  the  mismatch  5y/y»  and  as  a 
consequence  the  drift  of  the  particle  orbits  is,  for  all 
practical  purposes,  zero. 

The  electric  field  lines  and  field  components,  given  by  Eq. 
(2.33)  through  (2.37)  are  plotted  in  figs.  6b,  6c  and  6d.  The 
various  parameters  for  these  plots  are  listed  in  Table  III.  The 
dotted  lines  in  Fig.  6b  indicate  the  location  of  the  accelerating 
gap.  Fifty  one  terms  have  been  retained  in  the  series.  In  Fig. 
6c,  the  Gibbs  phenomenon  is  quite  noticeable 

at  p  =  0  =  10  cm,  where  the  E  component  is  discontinuous.  As 

W  ^ 

expected,  the  line  integral  along  the  s-axis  and  at  any 

radius  o  is  equal  to  V^.  As  may  be  seen  from  Fig.  6d 

the  E^  component  changes  sign  as  s  increases,  and  therefore,  on 

the  average  it  has  little  effect  on  the  particle  orbits.  The 

poloidal  component  E^  was  not  plotted  since  it  is  zero  at  4  »  o. 

In  general  it  is  small  and  its  contribution  to  the  particle 

orbi ts  is  negl i gi bl e . 

The  electric  fields  given  above  are  a  good  representation  of 
the  fields  produced  inside  the  torus  by  a  transmission  line, 
since  in  this  region  the  inductive  magnetic  field  is  zero  and 
therefore  the  potential  is  described  by  A  V  =»  0. 

III.  Simple  Theoretical'Model 

a.  Approximate  Equations  of  Motion  and  Conservation  Law 

The  exact  equations  of  motion  are  too  complicated  to  provide 
any  insight  into  the  effect  that  the  torsatron  and  vertical 
magnetic  fields  have  on  the  particle  orbits.  A  better 


understanding  is  gained  by  a  simple  model  which  provides  fairly 
accurate  results  and  is  based  on  the  assumption  that  the  toroidal 
velocity  v  of  the  particle,  being  very  close  to  the  velocity  of 
light,  remains  constant  (but  not  the  energy).  In  this  case 

s  =  Vgt  ,  (3.1) 

and  the  remaining  equations  of  the  motion  are 

(Vg  -  z  B^)  ,  (3.2a) 

(r  -  v^  B^).  (3.2b) 

For  simplicity,  the  self  fields  will  be  neglected  in  the  present 
analysis.  Also,  the  magnetic  field  of  the  torsatron  and  toroidal 
coils  is  given  by  the  simplified  Eqs.  (2.13)  and  that  of  the 
plates  by  Eqs.  (2.24).  In  Eqs.  (2.13),  2os  is  replaced 
by  -  (cf .  Eq.  (3.1) ), 

where 
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Set  X  *  r  -  r^,  and  define  the  complex  amplitude 


u  =  x  +  iz=oe  . 


1 


(3.4) 


Then  Eqs.  (3.2a),  (3.2b)  can  be  combined  into  one  complex 
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Equation  (3.5)  is  still  too  difficult  to  handle.  But  it 
leads  to  a  conservation  law  when  the  fast  motion  which  oscillates 
with  frequency  u  is  averaged  out  in  Eq.  (3.5).  Details  may  be 
found  elsewhere'  K  Only  the  final  result  will  be  given  here 

which  was  obtained  under  the  assumptions  that  y  remains  constant 
and  that  <<  |  Under  these  conditions  the 

conservation  law  becomes 
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Here,  K  is  the  constant  of  integration  and  is  determined  from  the 
initial  conditions,  (X,Z)  is  the  average  (or^^^iding  center) 


position  of  the  beam  centroid,  R  *  (X*"  +  1“)  ,  y^  is  the  value 

of  Y  that  matches  the  vertical  magnetic  field  from  the  plates 
when  the  beam  centroid  i s  at  X  *  Z  =  0  (i.e.,  at  the  origin  of 
the  local  coordinate  system),  and  5y/y  is  the  mismatch  between 


(X^  +  Z^) 


the  actual  value  of  y  and  y^.  When  the  mismatch  is  equal  to 
zero,  then  the  beam  is  matched  at  the  origin,  i.e.,  the  point  X  » 
Z  *  0  is  an  equilibrium  position. 

Consider  the  special  case  when  2oR  <<  1  and  Sy/y  <<  1. 

Then  Eq.  (3.8),  which  is  actually  an  equation  for  the  particle 
orbits,  simplifies  to  the  following  expression: 


(b  +  1  - 
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where  K  is  a  new  integration  constant.  This  equation  describes 
either  an  ellipse  or  a  hyperbola  centered  at  (Ar,0),  where 


AT.  a 

r  b  +  1  -  n.  ‘ 
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(3.11) 


At  the  initial  stage  of  the  acceleration,  the  torsatron 
magnetic  field  is  much  larger  than  the  vertical  magnetic  field 
B_-,  i.e.,  |B  /B,  1  >>  1.  Since  the  external  field  index  n.  is 

^  u  (  0  Z  0  *  ^ 

of  order  one,  it  follows  that  |b_j  is  much  greater  than  one,  or 
l^ol*  ’  ”ol’  simplifies  to  the  relation 


,  1^'|2  +  ,  1^'  ^  (3.12) 

'^0  '^o 

It  is  seen  then  that,  in  the  presence  of  the  torsatron  field,  at 
the  initial  stage  of  the  acceleration  the  orbits  on  the  (r,z) 
plane  are  circles  centered  at 


This  orbit  displacement  is  independent  of  the  external  field 
index  n^  and  since  |b  |  >>  1,  it  is  not  sensitive  to  the  value  of 
the  mismatch  5y/y,  i.e.,  the  torsatron  field,  which  dominates  at 
the  initial  stage  of  the  acceleration,  does  provide  confinement 
of  the  beam  inside  the  torus. 
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At  a  later  stage  of  the  acceleration,  as  the  vertical 
magnetic  field  increases  in  synchronism  with  the  energy  of 
the  beam,  it  becomes  eventually  greater  than  and  the 
condition  jb_[  <<  1  holds.  In  this  case,  Eq.  (3.10)  reduces  to 
the  expression 


-  %>  'h  -  *  "o 


(3.14) 


Therefore,  in  order  to  have  closed  orbits  during  the  whole 
duration  of  the  acceleration,  the  condition  0  <  n^<  1  must  be 
satisfied,  and  it  is  this  condition  that  requires  the  presence  of 
a  gap  on  the  outer  plate,  as  indicated  in  Section  II.  Also,  in 
this  case 


(3.15) 


which  is  sensitive  to  the  mismatch  fiy/y  and  to  the  external  field 
index  n^.  In  the  special  case,  when  n^^  *  0,  the  orbits  are  given 
by  the  relation  X  »  constant,  i.e.,  they  describe  a  vertical 
motion  and  this  may  provide  a  possible  means  of  extracting  the 
beam  from  the  toroidal  chamber. 

b.  Rate  of  Change  of  the  Current  in  the  Plates 

As  mentioned  in  the  previous  Sections,  the  two  coaxial 
cylindrical  plates  generate  a  local  vertical  magnetic  field  which 
increases  in  synchronism  with  the  beam  energy  so  that  the  major 
radius  of  the  ring  remains  approximately  constant.  From  Eq. 
(3.7),  a  particle  at  the  origin  of  the  local  coord inate 
system  e^,  e^,  e^,  will  gain  energy  at  the  rate 


e  f _ g_ 

me  2  irr 


^2  no 
c 


(3.16) 


Equation  (3.16)  has  been  derived  from  Eqs.  (2.15),  (2.26)  and  by 
taking  the  average  accelerating  field  in  the  gap  equal  to 
Egap  “  quantity  is  equal  to 
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The  beam  will  remain  matched  at  the  origin  during  the 
acceleration  when  (cf.  Eqs.  (3.9b),  (2.25b),  (2.25c)) 
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From  Eqs.  (3.16),  (3.18),  it  follows  that  the  rate  of  change  of 
the  current  in  the  plates  should  be  equal  to 
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The  parameter  np  is  the  acceleration  correction  factor  due  to  the 
elect ric  field  of  the  plates  and  has  been  listed  in  Table  II. 

It  is  interesting  to  show  that  even  when  the  beam  centroid 
is  slightly  off  the  origin  at  some  later  time,  the  mismatch  iSy  is 
still  zero  to  first  order  in  the  displacement  from  the  origin, 
if  Sy  was  initially  zero.  A  qualitative  argument  will  be 
presented  here  and  a  more  precise  proof  is  given  in  the 
Appendix.  If  the  beam  is  at  some  position  r  =  r^  +  x,  the 
average  accelerating  field  is  equal  to  ^ggp*  V^/(2iTr)  and  the 
vector  potential  at  that  position  is  proportional  to 

A  A 

{3A^/3x)  X,  where  the  partial  derivative  is  evaluated  at 

“0  n 

the  origin.  Therefore,  from  Eqs.  (3.9a),  (3.9b)  and  (3.7),  it 
follows  that 


J-  f 


0  t  irr , 


(3.20) 


In  deriving  Eq.  (3.20)  we  used  Eqs.  (3.19a)  and  (2.20b). 

It  is  seen  then  that,  in  the  rate  of  change  of  Sy,  associated 
with  the  slow  motion  of  the  beam,  the  first  order  effect  of  the 
plates  is  canceled  by  the  first  order  toroidal  corrections  of  the 
accelerating  field.  Therefore,  the  toroidal  corrections 
(especially  the  1/r  dependence)  of  the  accelerating  field  must  be 
incorporated  in  the  equations  of  motion.  Otherwise,  there  is  a 
slow  drift  of  the  beam  in  the  radial  direction.  This  has  been 
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observed  in  several  computer  runs  which  have  not  included  the 
toroidal  corrections.  After  the  toroidal  corrections  were 
included,  the  rate  of  change  of  5y  became  negligible. 


1  c.  Resonances  and  Multiple  Frequencies  in  the  Rebatron 

I 


In  contrast  to  the  modified  betatron  where  the  acceleration 
takes  place  continuously,  in  the  rebatron  it  is  a  stepwise 
process,  since  the  accelerating  field  is  localized  at  the  gap. 
Therefore,  there  is  a  periodicity  associated  with  the  rise  of  the 
beam  energy,  with  period  T  =  i.e.,  equal  to  the  time  it 

takes  the  particle  to  make  a  full  turn  around  the  torus.  Let  the 
increment  per  turn  in  y  be  Ay.  Since  the  increase  of  y  takes 
place  during  the  short  time  At  that  the  particle  crosses  the 
accelerating  gap.  At  <<  T,  and  y  is  a  ladder  function,  of  the 
step  Ay  and  length  T.  Therefore,  y  can  be  expressed  as  follows; 


Y(t) 


1  i  n 


(t)  + 


''osc"’’ 


(3.21a) 


where 


Ylin(^)  "  Yo  Ay  (y  +  y),  0  <t  <  »,  (3.21b) 

Yosc^^^  =  Ay  (^  -  f),  0  <  t  <  T,  (3.21c) 


^nd  Y^,,(t  +  T)  a  y^  (t).  Since  the  first  crossing  through 
os  c  os  c 

the  accelerating  gap  occurs  at  time  t  »  0  ,  after  m  turns  y  is 
equal  to  y  =  y^  +  (m  +  1)  Ay,  where  y  is  the  initial  value  ofy. 

Since  YQ5j.(t)  is  a  periodic  function,  it  can  be  Fourier 
analyzed ,  i.e.. 


A  similar  expression  has  been  obtained  for  a  periodic  disturbance 
in  the  modified  betat ron ^ ^ ,  The  difference  here  is  that  the 
occurrence  of  resonances  in  the  rebatron  is  an  inherent  propery 
of  the  device  due  to  the  localized  accelerating  electric  field 
and  the  periodic  nature  of  the  acceleration.  According  to  Eq. 
(3.9b),  resonance  occurs  each  time  y  ■  »  where 


^n  * 


2n^-  1 


(3.26) 


and  n  *  1,  2,  ...  As  it  will  be  shown  in  Section  IV,  these 
resonances  cause  the  orbits  to  expand  and,  in  the  worst  case 
(i.e.,  when  n  =  1)  the  beam  hits  the  wall  of  the  torus. 


It  is  apparent  from  Eq.  (3.24a)  that  always  (i)_  <  But, 


as  the  vertical  magnetic  field  increases, 
multiple  of  (o_,  i.e.,  *  mtii_  .  In  ten 
when 


may  become  a 


1  _  °zo 

FT 


m  -  1 


In  terms  of  e|  this  happens 


(3.27) 


where  m  =  2,  3,  ...  The  occurence  of  a  multiple  frequency 
changes  the  orbits  into  triangular  or  square  shapes  depending  on 
the  value  of  m.  According  to  Eq.  (2.9b)  a  multiple  frequency 


occurs  each  time  y  »  y^,  where 


I'^eol 


m  -  1 


(3.28) 


and  m  =  2,  3,  ...  This  has  been  observed  in  the  computer  runs 
and  examples  will  be  given  in  Section  IV. 


IV.  Numerical  Results  and  their  Interpretation 


The  various  properties  of  the  rebatron  were  investigated  by 
integrating  the  relativistic  equations  of  motion  for  the  beam 
centroid  using  Eqs.  (2.1)  through  (2.8)  for  the  torsatron  and 
external  toroidal  magnetic  fields,  Eqs.  (2.19)  through  (2.23)  for 
the  magnetic  field  of  the  plates,  Eq.  (2.26)  for  the  electric 
field  of  the  plates,  and  Eqs.  (2.33)  through  (2.37)  for  the 
accelerating  electric  field.  The  various  parameters  and  their 
values  are  listed  in  Tables  I,  II,  and  III.  The  self  fields  of 
the  beam  have  not  been  included  in  the  computations. 

To  test  of  the  extent  over  which  the  external  field  index 
confines  effectively  the  particle  orbits,  some  computer  runs 
without  acceleration  at  y  *  387  were  made  in  the  presence  of  a 
gap  on  the  outer  plate.  The  vertical  magnetic  field  was  set 
equal  to  »  6597  Gauss  (i.e.,  it  matches  the  value 
of  Yg  “  387).  The  results  for  various  initial  conditions  are 
shown  in  Fig.  7.  As  expected,  when  the  centroid  is  initially  at 
the  origin,  it  remains  there,  since  y  and  B,  are  matched  at  that 
position.  On  the  other  hand,  since  the  external  field  index 
depends  on  the  position  and  is  close  to  one  at 
X  »  6  cm,  z  *  0  (cf.  Fig.  5a),  when  the  particle  is  placed 
initially  there,  it  moves  in  a  radial  outward  direction  (cf.  Eq. 
(3.14))  and  eventually  strikes  the  wall.  Therefore,  for  the 
parameters  chosen,  the  particle  is  successfully  confined  Inside 
the  torus  up  to  distances  of  5  cm  from  the  minor  axis. 

In  order  to  test  the  focusing  effect  of  the  torsatron 
magnetic  field  at  the  initial  stage  of  the  acceleration  a  set  of 
runs  was  made  without  a  gap  on  the  outer  plate  so  that 
n^m  -  0.18.  The  results  are  shown  in  Fig.  8.  The  initial 
vertical  magnetic  fields  are:  a)  B^^  =  341  Gauss  (it 
matches  Yq  *  20)  ,  in  Fig.  8a,  b)  B^^  =  3.41  kG  (it  matches 
Yq  »  200)  in  Fig  8b,  and  c)  8^^  *  17.05  kG  (it 
matches  Yq"  1000)  in  Fig.  8c.  In  all  three  cases  the  rate  of 
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I 


change  of  the  current  in  the  plates  was  equal 

to  Ip  “  1.79  X  lO^^Amp  /sec.  The  linearly  increasing  vertical 
magnetic  field  of  the  plates  matches  the  dotted  lines  of  y  versus 
time  in  Figs.  8a",  8b"  and  8c".  It  is  seen  then  that  at  low 
values  of  y  the  torsatron  field  confines  the  orbit  inside  the 
torus  even  though  the  external  field  index  has  a  negative 
value.  On  the  other  hand,  at  high  values  of  y  the  torsatron 
field  is  no  longer  effective  in  confining  the  particle  orbit 
which  becomes  sensitive  to  the  external  field  index  (cf.  Eqs. 
(3.14),  (3.15))  . 

Figures  9  and  10  provide  the  main  results  of  this  paper. 

They  correspond  to  two  different  initial  positions  of  the 
particle  orbit.  The  various  parameters  for  these  computer  runs 
are  listed  in  Table  IV.  In  both  cases  the  outer  plate  had  a  gap 
so  that  n^a  0.5  at  the  minor  axis.  Figure  9A  shows  magnified  the 
particle  orbit  during  the  initial  steps  of  the  acceleration.  The 

dotted  line  in  Fig.  9A"  shows  the  value  of  y  versus  time  which 

matches  the  linearly  increasing  vertical  magnetic  field  in  the 
plates.  In  both  runs,  a  value  of  y  in  excess  of  1850  was 
reached  within  a  time  of  10  usec,  or  within  500  revolutions.  In 
both  runs  there  is  no  appreciable  drift  of  the  particle  orbits. 

As  stated  in  Section  III  b,  this  is  due  to  the  fact  that  the  rate 

of  change  of  the  mismatch  6y  has  a  second  order  dependence  on  the 
displacement  and  therefore,  if  4y  is  initially  zero,  it  remains 
vary  small  during  the  acceleration.  In  the  initial  computer 
runs,  the  toroidal  corrections  of  the  accelerating  field  were  not 
included  in  the  computations  and  a  very  slow  inward  or  outward 
radial  drift  as  well  as  an  extreme  sensitivity  to  the  initial 
value  of  the  vertical  magnetic  field  and  its  rate  of  change  were 
observed.  Values  of  y  beyond  1000  could  not  be  reached  since  the 
orbits  would  slowly  drift  and  hit  the  wall  of  the  torus.  After 
the  toroidal  corrections  were  included  in  the  computations,  the 
drift  disappeared  as  Figs.  9  and  10  indicate. 
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Another  feature  that  is  observed  in  Figs.  9  and  10  is  the 
expansion  of  the  orbits  as  time  evolves.  This  is  not  a 
continuous  expansion  but  it  occurs  rather  in  bursts,  the  most 
pronounced  being  when  y  exceeds  the  value  of  1300.  As  discussed 
in  Section  III  c,  in  the  rebatron  there  can  be  resonances  (cf.Eq. 
(3.26))  which  for  the  parameters  chosen  here,  occur 
at  Yj  *  1830,  y2  *  523,  Y3*  323,  y^  =  236,  etc.  The  resonance  at 
Y  *  1830  can  be  seen  clearly  in  Figs.  9f  and  lOj  and  causes 
eventually  the  orbit  to  hit  the  wall  of  the  torus.  This  is  the 
worst  resonance  and  occurs  when  l®9ol^®zo  *  The  amplitude 

of  the  nth  resonance  is  proportional  to  1/n  (cf.  Eq.  (3.23)), 
i.e.,  it  decreases  rather  slowly  with  n.  Therefore,  the  higher 
resonances  are  also  noticeable.  The  second  resonance  may  be  seen 
in  Fig.  11,  which  is  the  same  as  Fig.  9c,  but  magnified. 
Initially,  when  y  varies  from  500  to  600  the  orbit  rapidly 
expands  from  0.7  cm  to  2.2  cm,  but  as  y  keeps  increasing  from  600 
to  900  the  radial  position  of  the  orbit  remains  close  to  2.5  cm, 
i.e.,  it  does  not  vary  as  much.  A  comparison  of  the  size  of  the 
orbits  in  Figs.  10b,  10c  and  lOd,  indicates  that  there  is  a 
resonance  close  to  y  TT  500  in  the  second  run  too.  The  higher 
resonances  beyond  the  second  are  clustered  much  closer  together, 
their  amplitudes  are  smaller  and,  therefore,  they  are  not  as 
noticeable  individually.  But  they  do  cause  an  expansion  of  the 
particle  orbit  as  indicated  in  Fig.  9a  and  in  Figs.  10a  and 
10b.  The  value  of  the  total  toroidal  magnetic  field  used  in 
these  runs  was  ®  -  15.6  k6.  A  higher  value  of  B„  would 

cause  the  first  resonance  to  occur  at  a  value  of  y  greater  than 
1830  (cf.  Eq.  (3.26)).  Therefore,  there  is  no  intrinsic 
difficulty  in  the  rebatron  to  accelerate  the  beam  to  one  GeV  and 
beyond . 

Figure  9f  indicates  a  radical  change  in  the  shape  of  the 
particle  orbit  as  y  surpasses  the  value  of  1800.  A  similar 
observation  holds  for  Fig.  lOj  and  also  for  Fig.  lOf,  where  y  is 
in  the  neighborhood  of  1100.  As  stated  in  Section  IIIc,  this  is 
due  to  the  fact  that  the  characteristic  frequencies  become 
multiples  of  each  other.  For  the  parameters  chosen  here. 


Eq.  (3.28)  predicts  that  multiple  frequencies  should  occur 
when  y2  "  1830,  =  1121,  etc.  Notice  that  in  the  first  run 

there  is  no  change  at  y  =  1830.  This  indicates  that,  in 
addition  to  the  occurrence  of  a  multiple  frequency,  some  other 
conditions  must  be  satisfied  in  order  to  change  the  shape  of  the 
orbit.  A  clear  demonstration  is  given  in  the  runs  shown  in  Figs. 
12a  and  12b,  where  there  is  no  acceleration.  In  both  runs  the 
initial  position  is  =  2  cm,  =  0,  and  the  initial  velocity 
is  v^  =  v_  =  0,  v«  =  c.  Also,  in  both  runs,  n„  =  1/2.  In  the 

run  of  Fig.  12a,  y  =  1830  and  31.167  kG  (it  matches 

Y  =  1830),  while  in  the  run  of  Fig.  12b,  y  =  1121  and  B^g  = 
19.095  kG  (it  matches  y^  =  1121) 

V.  Conclusions 

A  detailed  numerical  and  analytical  investigation  of  the 
beam  dynamics  in  the  rebatron  accelerator  has  been  carried  out. 

It  has  been  shown  that  energies  approaching  one  GeV  can  be 
achieved  within  10  microseconds.  Since  the  acceleration  occurs 
in  such  a  short  time,  the  device  may  not  be  sensitive  to  the 
various  inabilities  and  the  loss  to  synchroton  radiation  should 
be  small. 

As  a  result  of  the  periodic  nature  of  the  acceleration  in 
the  rebatron,  it  has  been  found  that  resonances  can  be  excited 
which  cause  the  orbits  to  expand  and,  in  the  worst  case,  the  beam 
may  hit  the  wall  of  the  toroidal  chamber.  The  worst  resonance  is 
the  fundamental,  i.e.,  n  =  1,  and  it  occurs  when  the  vertical 
magnetic  field  becomes  twice  the  total  toroidal  magnetic  field. 
This  could  be  a  potential  problem  in  the  rebatron  unless  the 
toroidal  magnetic  field  is  sufficiently  high  so  that  the 
fundamental  resonance  is  not  excited.  In  the  cases  studied  here, 
the  total  toroidal  magnetic  field  was  set  at  15.6  KG  which 
allowed  the  beam  to  be  accelerated  at  a  y  slightly  higher  than 
1800  before  the  fundamental  resonance  was  reached 
at  Yj  “  1830.  From  the  computer  runs,  it  has  been  concluded 


that  the  higher  resonances  do  cause  the  orbits  to  expand  but  they 
are  not  as  detrimental  as  the  fundamental  resonance  is. 

To  keep  the  inductive  voltage  managable,  the  vertical  plates 
which  generate  the  vertical  magnetic  field  should  have  a  low 
inductance.  As  an  example,  for  the  parameters  in  table  II,  the 
inductance  Lp  of  the  plates  is  approxi mately 2.5  uH.  A 
vertical  magnetic  field  =  17  kG  is  required  to  reach  a  v  of 
1000,  while  the  plates  generate  at  the  minor  axis  17.58  Gauss  per 
kA.  Therefore,  the  current  in  each  plate  required  to  generate  17 
kG  is  Ip  »  967  kA.  Suppose  that  the  acceleration 
time  At  =  10  usec.  Then  the  rate  of  change  of  the  current  in 

the  plates  is  Ip  *  0.967  x  10^^  Amp/sec  and  the  inductive  voltage 
is  Vp  »  242  kV,  which  is  fairly  high. 

The  simple  analytical  model  of  Section  III,  provided  quite 
reliable  results  with  regard  to  the  confinement  of  the  beam  by 
the  torsatron  magnetic  field.  The  predictions  of  the  model  are 
in  good  agreement  with  the  numerical  results.  It  has  been  shown 
that,  at  the  initial  stage  of  the  acceleration,  the  torsatron 
field  provides  excellent  confinement  of  the  beam,  and  the  beam 
equilibrium  position  was  quite  insensitive  to  the 
mismatch  Sy/y  and  to  the  external  field  index  n^,  as  long 
as  0)^  I  beam  energy  increases,  the  stage  is 

reached,  where  I torsatron  field  is  no  longer 

effective  and  the  beam  equilibrium  position  becomes  sensitive  to 
both  the  mismatch  Sy/y  and  the  external  field  index  n^.  As 
mentioned  in  Section  III  b,  the  sensitivity  in  n^  may  prove 
beneficial  in  providing  a  means  for  extracting  the  beam.  (cf. 

F  i  gs .  8b  and  8c ) . 

Another  interesting  feature  of  the  acceleration  in  the 
rebatron  is  the  absence  of  any  drift  of  the  beam  when  the 
vertical  magnetic  field  is  matched  to  the  beam  energy.  As 
explained  in  Section  III  b,  this  is  due  to  the  fact  that,  the 
first  order  toroidal  correction  of  the  accelerating  field  cancels 
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out  the  first  order  correction  of  the  electric  field  from  the 
plates.  Therefore,  if  the  mismatch  is  initially  zero,  it  will 
remain  very  small  during  the  acceleration,  and  there  is  no 
drift.  This  has  been  verified  in  several  computer  runs. 

In  the  analytical  and  numerical  results  presented  in  this 
paper,  the  self  fields  of  the  beam  have  not  been  included.  Thei 
effect  on  the  beam  dynamics  and  the  acceleration  process  will  be 
the  subject  of  a  future  publication. 


Table  II 
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Parameters  associated  with  the  magnetic  field  of  the  plates  shown 
In  Figs.  4  and  5. 


Torus  major  radius  r  (cm) 


Toroidal  chamber  minor  radius  p,^(cm) 


Inner  plate  distance  from  minor  axis  a^  (cm) 
Half  height  of  outer  plate  hg/2  (cm) 

Half  height  of  outer  gap  Hg/2  (cm) 

Half  height  of  Inner  plate  h^/2  (cm) 


Current  In  outer  plate  Ip  (kA) 


at  minor  axis  (Gauss) 


External  field  Index  at  minor  axis  n^ 
Acceleration  correction  factor 


Fig  4 

Fig  5 

=  100 

100 

*  10 

10 

(cm) 

=  12 

12 

(cm) 

*  -12 

-12 

»  24 

24 

»  0 

»  24 

*  1 


2.35 

24 

1 


«  18.45  17.58 

*  -0.18  0.5 

»  1.038  1.035 


Table  IV 

Parameters  of  the  runs  shown  1n  Figs.  9  and  10. 


Initial  current  in  outer  plate  Ip  (kA) 

Initial  vertical  magnetic  field  from 
plates  B^g  (Gauss) 

Initial  matching  value  of  relativistic 
factor 

m 

Rate  of  change  of  current  in  plates 
Ip  (Amp/sec) 

Accelerating  voltage  (MV) 

Increment  of  y  per  revolution 
Initial  relativistic  factor  y 
Initial  positions  (cm)  r=100 


=  8.1 


*  142.4 


»  8.35 


»  1.8745 


*  -2 


»  3.914 

»  7 


z=0  Fig.  9 
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z  =  1.0  Fig.  10 
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Initial  velocities 
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Fig.  3  —  Torsatron  magnetic  field  components  (a)  at  <I>  -•  s  ■  o  and  (b)  at  s  —  o,  d*  —  jr/2.  In 
addition  to  the  torsatron  field,  there  is  a  toroidal  field  — 10  kG  which  is  generated  by  a  set  of 

toroidal  coils,  (c)  Magnetic  field  lines  in  (r,z)  and  (p,s)  planes. 


(Gauttal 


p  (cm) 

Fig.  3  (Cont’d)  —  Torsatron  magnetic  field  components  (a)  at  1>  —  s  *  o  and  (b)  at  a  toroidal  field 
BT*  -  — 10  kG  which  is  generated  by  a  set  of  toroidal  coils,  (c)  Magnetic  field  lines  in  (r,z)  and  (p,s) 
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Fig.  4  —  Magnetic  field  of  plates  without  a  gap  on  the  outer  plate,  at  various  distances  along  the 
r-axis  or  the  z-axis.  (a)  B.  versus  r,  (b)  B,  versus  r,  (c)  B,  versus  z,  and  (d)  Bj  versus  z. 
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Effect  of  the  external  field  index  n^  of  the  magnetic  field  of  the  plates  on  the  particle  orbit 
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Fig.  8  —  Demonstration  that  the  particle  orbit  is  confined  by  the  torsatron  fields  only  at  the  initial  stage 
of  the  acceleration.  In  all  cases,  there  is  no  gap  on  the  outer  plate,  so  that  n,,  11  -  0.18.  Fig.  (7a)  is 
the  same  as  Fig.  (7a)  but  on  a  magnified  scale. 
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Fig.  3  (Cont’d)  —  Demonstration  that  the  panicle  orbit  is  confined  by  the  torsatron  fields  only  at  the 
initial  stage  of  the  acceleration.  In  all  cases,  there  is  no  gap  on  the  outer  plate,  so  that  ng  n  -  0.18. 
Fig.  (7 A)  is  the  same  as  Fig.  (7a)  but  on  a  magnified  scale. 
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Fig.  8  (Corn’d)  —  Demonstration  that  the  particle  orbit  is  confined  by  the  torsatron  fields  only  at  the 
initial  stage  of  the  acceleration.  In  all  cases,  there  is  no  gap  on  the  outer  plate,  so  that  n,  11  -  0.18. 
Fig.  (7A)  is  the  same  as  Fig.  (7a)  but  on  a  magnified  scale. 


GAMMA  z  (cm) 


Partici*  Orbit  in  the  Transverse  Plane 


r  (cm) 


1013 

1011 

1009 

1007 

1005 

1003 

1001 


0 


20  40 

TIME  (nsec) 


Fig.  8  (Cont’d)  —  Demonstration  that  the  particle  orbit  is  confined  by  the  torsatron  fields  only  at  the 
initial  stage  of  the  acceleration.  In  all  cases,  there  is  no  gap  on  the  outer  plate,  so  that  Og  11  —  0.18. 
Fig.  (7A)  is  the  same  as  Fig.  (7a)  but  on  a  magnified  scale. 
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Fig.  9  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  9(a)  through  9(f))  and  the  corresponding  relativistic 
factor  y  (Figs.  9 (a')  through  9(r)).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z  —  0  cm. 
The  particle  orbit  and  y  at  the  initial  stage  of  the  acceleration  are  shown  on  Figs.  9(A)  and  9(A')). 
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Fig.  9  (Cont’d)  —  Particle  orbit  on  the  (r,2)  plane  (Figs.  9(a)  through  9(f))  and  the  corresponding 
relativistic  factor  y  (Figs.  9(a')  through  9(f)).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z 
~  0  cm.  The  particle  orbit  and  y  at  the  initial  stage  of  the  acceleration  are  shown  on  Figs.  9(A)  and 
9(A')). 
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Fig.  9  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  9(a)  through  9(f))  and  the  corresponding 
relativistic  factor  y  (Figs.  9(a')  through  9(r)).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z 
*  0  cm.  The  particle  orbit  and  y  at  the  initial  stage  of  the  acceleration  are  shown  on  Figs.  9(A)  and 
9(A')). 
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Fig.  9  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  9(a)  through  9(f))  and  the  corresponding 
relativistic  factor  y  (Figs.  9(a')  through  9(r)).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z 
“  0  cm.  The  particle  orbit  and  y  at  the  initial  stage  of  the  acceleration  are  shown  on  Figs.  9(A)  and 
9(A')). 
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Fig.  9  (Cont’d)  —  Panicle  orbit  on  the  (r,z)  plane  (Figs.  9(a)  through  9(f))  and  the  corresponding 
relativistic  factor  y  (Figs.  9(a')  through  9(f)).  The  initial  position  of  the  particle  is  at  r  *  100  cm,  z 
“  0  cm.  The  particle  orbit  and  y  at  the  initial  stage  of  the  acceleration  are  shown  on  Figs.  9(A)  and 
9(A')). 
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Fig.  9  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  9(a)  through  9(f))  and  the  corresponding 
relativistic  factor  y  (Figs.  9(a')  through  9(r)).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z 
*■  0  cm.  The  particle  orbit  and  y  at  the  initial  stage  of  the  acceleration  are  shown  on  Figs.  9(A)  and 
9(A')). 
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Fig.  10  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding  relativistic 
factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  is  at  r  *  100  cm,  z  *  1  cm. 
The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(]))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  is  at  r  *  100  cm,  z 
"  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  is  at  r  “  100  cm,  z 
*  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  s  at  r  —  100  cm,  z 
—  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  is  at  r  -■  100  cm,  z 
■  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 


10 


Particle  Orbit  in  the  Transverse  Plane 


1160- 

1120- 

<  1080- 
2 
2 
< 
a 

1040- 

1000- 


960  - - ^ - 

5.0  5.2  5.4  5.6  5.8  6.0 

TIME  (^sec) 

Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  is  at  r  *■  100  cm,  z 
*  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(r)).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z 
*  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Panicle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  lOO")).  The  initial  position  of  the  particle  is  at  r  —  100  cm,  z 
“  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  lOl/)).  The  initial  position  of  the  particle  is  at  r  ~  100  cm,  z 
■■  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Fig.  10  (Cont’d)  —  Particle  orbit  on  the  (r,z)  plane  (Figs.  10(a)  through  10(j))  and  the  corresponding 
relativistic  factor  y  (Figs.  lO(a')  through  10(j')).  The  initial  position  of  the  particle  is  at  r  “  100  cm,  z 
»  1  cm.  The  parameters  of  this  run  are  listed  in  Table  IV. 
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Appendi  x 


In  the  modified  betat^on^^'  the  vector  potential  is  given  by 
the  exp ress i on ^ ^ 


A 


9 


(Al) 


where  n  is  the  external  field  index  and  r^  is  the  major  radius 
of  the  torus.  Set  r  *  r^  +  x  and  expand  Eq.  (Al)  close  to  the 
origin  x  *  z  =  0.  Then,  the  vector  potential  is  approximately 
equal  to 


S  "“zo  --o  (f  )  1. 

0  0 


and  the  rate  of  change  in  y  becomes: 


(A2) 


e  9 


me 


2  at 


r,  fl  .  J  (-^)S 

me  0  0 


(A3) 


The  same  expression  will  be  shown  to  be  true  in  the  rebatron  for 
the  slow  motion.  Notice  that,  because  of  the  toroidal  structure 
of  the  device,  the  accelerating  electric  field  is  periodic  with 
respect  to  the  variable  »  -  r^9  ,  with  period  2irr^. 

Therefore,  the  Eg  component  can  be  Fourier  analyzed,  i.e., 

ca 

E5(x,z,rg9)  »  A^(x,z)  +  2rA^(x,z)cos  (n9)  +  B  ^  ( x,  z  )  s  i  n  (n9 )  1 .  (A4) 

n  »  1 


where 


4*3 
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Ag(X,Z) 


It 

27  * 


(A5a  ) 


A^(x,z) 


TT 

7 J EgCx.z.r^e)  cos(n9)  de  , 


(A5b) 


B„(x.z) 


IT 

i/ 


^^(x.z.rQS)  sin(n9)  de. 


{A5c) 


and  9  *  -  0  .  The  position  at  which  the  electric  field 
component  is  computed  lies  at  a  fixed  distance  ^ 

from  the  major  axis.  Therefore,  in  order  to  change  the  Integrals 
over  an  angle  into  line  integrals,  the  transformation  9  *  s/r 
should  be  made.  In  particular,  under  this  t ransf ormat i on ,  Eq. 
(A5a)  becomes: 


irr 


(A6) 


Since  E^  is  localized  in  the  region  of  the  accelerating  gap,  the 
line  integral  above  is  equal  to  i.e.. 


Ag(x,z ) 


(A7) 


The  rate  of  change  in  y  is  approximately  equal  to 


f  E  +  ^ 


(A8) 


r,i 


m9, 


Av;' 


Replace  from  Eqs.  (A4),  (A7)  and  take  into  account 
that  ?  TT  -  (D^t,  where  2Tr/a)^  is  the  time  period  for  each 
revolution  around  the  torus.  Then  Eq.  (A8)  can  be  written  as 
f 0 1 1 ows : 


Y  =  Y 


where 


slow  *  ^osc  * 


(A9) 


s  1  ow 


e  f  _o 


me  '  2irr  c  5t 


.  i 


(AlOa) 


cosCnai^t)  -  B^(x,z)  sin(nti)jjt)'l .  (AlOb) 


The  slow  motion  or  the  drift  of  the  beam  is  associated  with  the 
first  term  in  Eq.  (A9)  while  the  resonances,  discussed  in  Section 
IIIc,  are  associated  with  the  second  term. 

Concentrate  on  the  slow  motion  and  use  Eqs.  (2.15),  (2.19a) 
and  (2,19b).  Then  Eq,  (AlOa)  becomes: 


slow  me 


and  also 


"  2??  ^  (x,z)'). 


(All) 


1  * 

T-^  +  -  A^ 
3x  re 


(A12a) 


(A12b) 


Take  a  Taylor  expansion  of  A^(x,z)  at  x=0,  z=0,  i.e.. 


3A  3A  ,  3^A  3^A  2  ^ 

%(x.z)  n  _  X  t  ^  z  t  ^  ^  ^  z  ?  xz, 

3  X  3  Z 


(A13) 


where  the  partial  derivatives  are  computed  at  x»0,  z*0.  It 
follows  from  Eqs.  (2.23),  (2.22a),  (2.22b)  and  (2.18)  that 
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3x  3x3z 


0  . 


(A14) 


Also  use  the  identities: 
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(A15b) 


Then  Eq.  (All)  becomes  approximately  equal  to 
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S»;Sv 


where  1/r  =  expanded  in  powers  of  x.  For  a 

matched  beam,  Ip  is  given  by  Eqs.  (3.19a),  (3.19b)  which  can  also 
be  written  as  follows  (cf.  Eq.  (A12b)): 


(A17) 


Therefore,  Eq.  (A16)  simplifies  to  the  expression 


^low-  -T  ®zo^^^  *  7  *  J-  -  7  T~  Jz- 

me  0  0 


(A18) 


If  the  linearized  expressions  for  B^,  B^  are  used  (cf.  Eqs 
(2.24a),  (2.24b)),  then  Eq.  (A18)  becomes: 


slow 


1  *  n 


0  0 


(A19) 


which  is  identical  to  Eq.  (A3)  for  the  modified  betatron. 
Therefore,  the  rate  of  change  of  the  mismatch  5y,  associated  with 
the  slow  motion  in  the  rebatron,  is  a  second  order  effect  in  the 
displacement  of  the  beam,  that  is  very  small  and  if  the 
mismatch  -Sy  is  initially  zero  is  will  remain  negligible. 
Coneequent ly ,  the  drift  of  the  beam  is  negligible  as  observed  in 
the  computer  runs. 
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